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O . Abstract 

Q^ I First inclusive measurements of isolated prompt photons in photoproduction at 

^ ' the HERA ep collider have been made with the ZEUS detector, using an integrated 

^ . luminosity of 38.4 pb~^. Cross sections are given as a function of the pseudorapidity 

Q^l and the transverse energy (77'*', Ej.) of the photon, for Ej^ > 5 GeV in the jp centre- 

of-mass energy range 134-285 GeV. Comparisons are made with predictions from 
Monte Carlo models having leading-logarithm parton showers, and with next-to- 
leading-order QCD calculations, using currently available parameterisations of the 
p\ ' photon structure. For forward r]^ (proton direction) good agreement is found, but 

j^ . in the rear direction all predictions fall below the data. 
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1 Introduction 

One of the primary aims of photoproduction measurements in ep collisions at HERA is 
the elucidation of the hadronic behaviour of the photon. The measurement of jets at 
high transverse energy has provided much information in this area 0, @. In the study 
of inclusive jets, next-to-leading order (NLO) QCD calculations are able to describe the 
experimental data over a wide range of kinematic conditions, although the agreement 
is dependent on the jet algorithm 0. However, significant discrepancies between data 
and NLO theories are found in dijet measurements [Q. A further means to study pho- 
toproduction is provided by final states with an isolated high-transverse-energy photon. 
These have the particular merit that the photon may emerge directly from the hard QCD 
subprocess ("prompt" photons), and also can be investigated without the hadronisation 
corrections needed in the case of quarks or gluons. In a previous measurement by ZEUS 
at HERA [^, it was shown that prompt photons, accompanied by balancing jets, are 
produced at the expected level in photoproduction and with the expected event charac- 
teristics. This work is extended in the present paper through the use of a much larger 
event sample taken in 1996-97, corresponding to an integrated ep luminosity of 38.4 pb~^. 
This allows a measurement of inclusive prompt photon distributions as a function of pseu- 
dorapidity rj'^ and transverse energy Ej, of the photon, and a comparison with LO and 
NLO QCD predictions. 



2 Apparatus and trigger 

During 1996-97, HERA collided positrons with energy Ee = 27.5 GeV with protons of 
energy Ep = 820 GeV. The luminosity was measured by means of the bremsstrahlung 
process ep — > e'jp. 

A description of the ZEUS apparatus and luminosity monitor is given elsewhere [^ . Of 
particular importance in the present work are the uranium calorimeter (CAL) and the 
central tracking detector (CTD). 

The CAL ^ has an angular coverage of 99.7% of 47r and is divided into three parts (FCAL, 
BCAL, RCAL), covering the forward (proton direction), central and rear angular ranges, 
respectively. Each part consists of towers longitudinally subdivided into electromagnetic 
(EMC) and hadronic (HAC) cells. The electromagnetic section of the BCAL (BEMC) 
consists of cells of ~20 cm length azimuthally and mean width 5.45 cm in the Z direction^, 
at a mean radius of ~1.3 m from the beam line. These cells have a projective geometry as 
viewed from the interaction point. The profile of the electromagnetic signals observed in 
clusters of cells in the BEMC provides a partial discrimination between those originating 
from photons or positrons, and those originating from neutral meson decays. 

The CTD is a cylindrical drift chamber situated inside a superconducting solenoid 
which produces a 1.43 T field. Using the tracking information from the CTD, the vertex 
of an event can be reconstructed with a resolution of 0.4 cm in Z and 0.1 cm in X, Y. In 



^ The ZEUS coordinate system is right-handed with positive-Z in the proton beam direction and an 
upward-pointing Y axis. The nominal interaction point isa,tX = Y = Z = 0. 
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this analysis, the CTD tracks are used to reconstruct the event vertex, and also in the 
selection criteria for high-ET photons. 

The ZEUS detector uses a three-level trigger system, of which the first- and second-level 
triggers used in this analysis have been described previously |^. The third-level trigger 
made use of a standard ZEUS electron finding algorithm ^ to select events with an 
electromagnetic cluster of transverse energy E^ > 4 GeV in the BCAL, with no further 
tracking requirements at this stage. These events represent the basic sample of prompt 
photon event candidates. 



3 Event selection 

The offline analysis was based on previously developed methods 0. An algorithm for 
finding electromagnetic clusters was applied to the data, and events were retained for 
final analysis if a photon candidate with Et > 5 GeV was found in the BCAL. A photon 
candidate was rejected if a CTD track, as measured at the vertex, pointed to it within 0.3 
radians; this removed almost all high-E'^ positrons and electrons, including the majority 
of those that underwent hard radiation. The BCAL requirement restricts the photon 
candidates to the approximate pseudorapidityQ range —0.75 < r]'^ < 1.0. 

Events with an identified deep inelastic scattered (DIS) positron in addition to the BCAL 
photon candidate were removed, thus restricting the acceptance to incident photons of 
virtuality Q^ < 1 GeV^. The quantity y™-^"'' = J2{.E — pz)/2Ee was calculated, where 
the sum is over all calorimeter cells, E is the energy deposited in the cell, and pz = 
Ecos9. When the outgoing positron is not detected in the CAL, y"^'^"-'^ is a measure 
oi y = Ey^in/Ee, where E^^in is the energy of the incident photon. If the outgoing 
positron is detected in the CAL, y"^'^"-^ ^ 1. A requirement of 0.15 < y"^^"-"^ < 0.7 was 
imposed; the lower cut removed some residual proton-gas backgrounds while the upper cut 
removed remaining DIS events, including any with a photon candidate that was actually 
a misidentified DIS positron. Wide-angle Compton scattering events {ep -^ e7p) were 
also excluded by this cut. This range of accepted y'^'^"''^ values corresponds approximately 
to the true y range 0.2 < y < 0.9. 

An isolation cone was imposed around the photon candidate: within a cone of unit radius 
in {r], (p), the total Et from other particles was required not to exceed O.lii^j^. This was 
calculated by summing the Et in each calorimeter cell within the isolation cone. Further 
contributions were included from charged tracks which originated within the isolation 
cone but curved out of it; the small number of tracks which curved into the isolation cone 
were ignored. The isolation condition much reduces the dijet background by removing 
a large majority of the events where the photon candidate is closely associated with a 
jet and is therefore either hadronic (e.g. a 7r°) or else a photon radiated within a jet. In 
particular, the isolation condition removes most dijet events in which a photon is radiated 
from a final-state quark. Approximately 6000 events with Ej, > 5 GeV remained after 
the above cuts. 



^ All kinematic quantities are given in the laboratory frame. Pseudorapidity 77 is defined as 
— lntan(0/2), where 9 is the polar angle relative to the Z direction, measured from the Z position 
of the event vertex. 



Studies based on the single-particle Monte Carlo samples showed that the photon energy- 
measured in the CAL was on average less than the true value, owing to dead material in 
front of the CAL. To compensate for this, an energy correction, typically 0.2 GeV, was 
added. 



4 Monte Carlo simulations 

In describing the hard interaction of photons of low virtuality with protons, two major 
classes of diagram are important. In one of these the photon couples in a pointlike way 
to a qq pair, while in the other the photon interacts via an intermediate hadronic state, 
which provides quarks and gluons which then take part in the hard QCD subprocesses. 
At leading order (LO) in QCD, the pointlike and hadronic diagrams are distinct and are 
commonly referred to as direct and resolved processes, respectively. 

In the present analysis, three types of Monte Carlo samples were employed to simulate: 
(1) the LO QCD prompt photon processes, (2) dijet processes in which an outgoing quark 
radiated a hard photon (radiative events), and (3) single particles (7, 7r°, rf) at high Et- 
All generated events were passed through a full GEANT-based simulation ||T^ of the 
ZEUS detector. 



The PYTHIA 5.7 [TTI and HERWIG 5.9 [|l|] Monte Carlo generators were both used to 



simulate the direct and resolved prompt photon processes. These generators include LO 
QCD subprocesses and higher-order processes modelled by initial- and final-state parton 



showers. The parton density function (pdf) sets used were MRSA |jT3| for the proton, 
and GRV(LO) [|I^ for the photon. The minimum p^ of the hard scatter was set to 
2.5 GeV. No multi-parton interactions were implemented in the resolved samples. The 
radiative event samples were likewise produced using direct and resolved photoproduction 
generators within PYTHIA and HERWIG. 

In modelling the overall photoproduction process, the event samples produced for the 
separate direct, resolved and radiative processes were combined in proportion to their 
total cross sections as calculated by the generators. A major difference between PYTHIA 
and HERWIG is the smaller radiative contribution in the HERWIG model. 

Three Monte Carlo single-particle data sets were generated, comprising large samples of 
7, tt'^ and rj. The single particles were generated uniformly over the acceptance of the 
BCAL and with a flat E^ distribution between 3 and 20 GeV; ii^r-dependent exponential 
weighting functions were subsequently applied to reproduce the observed distributions. 
These samples were used in separating the signal from the background using shower 
shapes. 



5 Evaluation of the photon signal 

Signals in the BEMC that do not arise from charged particles are predominantly due to 
photons, 7r° mesons and tj mesons. A large fraction of these mesons decay into multiphoton 



final states, the vr'' through its 27 channel and the rj through its 27 and Stt^ channels. 
For 7r°, 7] produced with Et greater than a few GeV, the photons from the decays are 
separated in the BEMC by distances comparable to the BEMC cell width in Z. Therefore 
the discrimination between photons and neutral mesons was performed on the basis of 
cluster-shape characteristics, thus avoiding any need to rely on theoretical modelling of 
the background. 

A typical high-E'r photon candidate consists of signals from a cluster of 4-5 BEMC cells. 
Two shape-dependent quantities were used to distinguish 7, 7r° and rj signals 0. These 
were (i) the mean width <6Z> of the cell cluster in Z, which is the direction of finer seg- 
mentation of the BEMC, and (ii) the fraction fmax of the cluster energy found in the most 
energetic cell in the cluster. The quantity <6Z> is defined as Y. \Eceii\Zceii — Z\] /J2 Eceiu 
summing over the cells in the cluster, where Z is the energy- weighted mean Z value of the 
cells. The <5Z> distribution for the event sample is shown in Figure 1(a), in which peaks 
due to the photon and 7r° contributions are clearly visible.^ The Monte Carlo samples of 
single 7, 7r° and rj were used to establish a cut on <5Z> at 0.65 BEMC cell widths, such 
as to remove most of the 77 mesons but few of the photons and vr'^s. Candidates with lower 
<5Z> were retained, thus providing a sample that consisted of photons, 7r° mesons and 
a small admixture of r] mesons. 

The extraction of the photon signal from the mixture of photons and a neutral meson 
background was done by means of the fmax distributions. Figure 1(b) shows the shape 
of the fmax distribution for the final event sample, after the <5Z> cut, fitted to the r] 
component determined from the <6Z> distribution plus freely-varying 7 and 7r° contribu- 
tions. Above an fmax value of 0.75, the distribution is dominated by the photons; below 
this value it consists mainly of meson background. Since the shape of the fmax distribu- 
tion is similar for the 77 and 7r° contributions, the background subtraction is insensitive to 
uncertainties in the fitted ti^ to t] ratio. 

The numbers of candidates with fmax > 0.75 and fmax < 0.75 were calculated for the 
sample of events occurring in each bin of a measured quantity. From these numbers, and 
the ratios of the corresponding numbers for the fmax distributions of the single particle 
samples, the number of photon events in the given bin was evaluated S . 



6 Cross section calculation and systematic uncertain- 
ties 

Cross sections are given for the photoproduction process ep -^ 7(prompt) + X, taking 
place in the incident jp centre-of-mass energy {W) range 134-285 GeV, i.e. 0.2 < y < 0.9. 
The virtuality of the incident photon is restricted to the range Q^ <1 GeV^, with a 
median value of approximately 10"'^ GeV^. The cross sections represent numbers of 
events within a given bin, divided by the bin width and integrated luminosity. They 
are given at the hadron level, with an isolation cone defined around the prompt photon 



•^The displacement of the photon peak from the Monte Carlo prediction does not affect the present 
analysis; the poor fit in the region <SZ> — 0.6-1.0 is taken into account in the systematic errors. 



as at the detector level. To obtain the hadron-level cross sections, bin-by-bin correction 
factors were applied to the corresponding detector-level distributions; these factors were 
calculated using PYTHIA. 

The following sources of systematic error were taken into account: 

1. Calorimeter simulation: the uncertainty of the simulation of the calorimeter re- 
sponse [0 gives rise to an uncertainty on the cross sections of ±7%; 

2. Modelling of the shower shape: uncertainties on the agreement of the simulated fmax 
distributions with the data correspond to a systematic error averaging ±8% on the 
final cross sections; 

3. Kinematic cuts: the cuts defining the accepted kinematic range at the detector level 
were varied by amounts corresponding to the resolution on the variables. Changes 
of up to 5% in the cross section were observed; 

4. ri/{ri + vr'') ratio: the fitted value was typically 25%; variations of this ratio in the 
range 15-35% led to cross section variations of around ±2%; 

5. Vertex cuts: narrowing the vertex cuts to (—25, +15) cm from their standard values 
of (—50, +40) cm gave changes in the cross sections of typically ±4%. 

In addition, studies were made of the effects of using HERWIG instead of PYTHIA for the 
correction factors, of varying the Et distribution applied to the single-particle samples, 
and of varying the composition of the Monte Carlo simulation in terms of direct, resolved 
and radiative processes. These gave changes in the cross sections at the 1% level. The 
1.6% uncertainty on the integrated luminosity was neglected. The individual contributions 
were combined in quadrature to give the total systematic error. 



7 Theoretical calculations 

In presenting cross sections, comparison is made with two types of theoretical calculation, 
in which the pdf sets taken for both the photon and proton can be varied, although there 
is little sensitivity to the choice of proton pdf. These are: 

(i) PYTHIA and HERWIG calculations evaluated at the final-state hadron level, as out- 
lined in Sect. 4. Each of these programs comprises a set of LO matrix elements augmented 
by parton showers in the initial and final states together with hadronisation; 



(a) NLO parton-level calculations of Gordon (LG) |T^,|T^ and of Krawczyk and Zem 



brzuski (K&Z) [|T8|. Pointlike and hadronic diagrams at the Born level are included, 
together with virtual (loop) corrections and terms taking into account three-body final 
states. The radiative terms are evaluated by means of fragmentation functions obtained 
from experiment. In both calculations, the isolation criterion was applied at the parton 
level. 



The LG and K&Z calculations differ in several respects [p!q-[T9|]. The K&Z calculation 
includes a box-diagram contribution for the process 'jg — >• •yg |^, but excludes higher- 
order corrections to the resolved terms which are present in LG. A value of Ajj^ = 200 



MeV (5 flavours) is used in LG while in K&Z a value of 320 MeV (4 flavours) is used, so 
as to reproduce a fixed value of a^ = 0.118 at the Z^ mass. The standard versions of both 
calculations use a QCD scale of p^. Both calculations use higher-order (HO) versions of 
the GRV ig and GS pi photon pdf sets. 



8 Results 

Figure 2 and Table 1 give the inclusive cross-section da/dEj, for the production of isolated 
prompt photons in the range —0.7 < 7]'^ < 0.9 for 0.2 < y < 0.9. All the theoretical models 
describe the shape of the data well; however the predictions of PYTHIA and especially 
HERWIG are too low in magnitude. The LG and K&Z calculations give better agreement 
with the data. 

Figure 3 and Table 2 give the inclusive cross-section da /drp for isolated prompt photons 
in the range 5 < Ej^ < 10 GeV for 0.2 < y < 0.9. Using the GRV pdf's in the photon, 
PYTHIA gives a good description of the data for forward pseudorapidities. The HERWIG 
distribution, while similar in shape to that of PYTHIA, is lower throughout; this is 
attributable chiefly to the lower value of the radiative contribution in HERWIG (see Sect. 
4). The LG and K&Z calculations using GRV are similar to each other and to PYTHIA. 
All the calculations lie below the data in the lower r]"' range. 

The effects were investigated of varying some of the parameters of the K&Z calculation 
relative to their standard values (NLO, 4 flavours, Ap-^ = 320 MeV, GRV photon pdf). 
Reducing the number of flavours used in the calculation to three (with A^g^ = 365 MeV) 
reduced the cross sections by 35-40% across the t]"' range, confirming the need to take 
charm into account. A LO calculation (with ^-j^ = 120 MeV and a NLO radiative con- 
tribution) was approximately 25% lower than the standard NLO calculation. Variations 
of the QCD scale between 0.25E^ and 4E^ gave cross-section variations of approximately 
±3%. 

Figure 3(b) illustrates the effects of varying the photon parton densities, comparing the 



results using GRV with those using GS. The ACFGP parton set |2^ gives results (not 
shown) similar to GRV. All NLO calculations describe the data well for rj'^ > 0.1, as does 
PYTHIA, but are low at more negative i]'^ values, where the curves using the GS parton 
densities give poorer agreement than those using GRV. 

As a check on the above results, the same cross sections were evaluated with the additional 
requirement that each event should contain a jet (see 0) with Et > 5 GeV in the 
pseudorapidity range (—1.5, 1.8). Both the measured and theoretical distributions were 
found to be of a similar shape to those in Fig. 3. 

The discrepancy between data and theory at negative r]'^ is found to be relatively strongest 
at low values oiy. Figure 4 shows the inclusive cross section da /drp as in Fig. 3, evaluated 
for the three y ranges 0.2-0.32, 0.32-0.5 and 0.5-0.9 by selecting the y'^'^"-^ ranges 0.15- 
0.25, 0.25-0.4 and 0.4-0.7 at the detector level. The numerical values are listed in Table 
2. In the lowest y range, both theory and data show a peaking at negative rp , but it is 
stronger in the data. The Monte Carlo calculations indicate that the peak occurs at more 



negative r]"' values as y increases, eventually leaving the measurement acceptance. In the 
highest y range (Fig. 4(c)), agreement is found between theory and data. The movement 
of the peak can be qualitatively understood by noting that for fixed values of Et and x^, 
where x^ is the fraction of the incident photon energy that contributes to the resolved 
QCD subprocesses, measurements at increasing y correspond on average to decreasing 
values of pseudorapidity. By varying the theoretical parameters, the discrepancy was 
found to correspond in the K&Z calculation to insufficient high x^ partons in the resolved 
photon. 



9 Summary and conclusions 

The photoproduction of isolated prompt photons within the kinematic range 0.2 < ?/ < 
0.9, equivalent to incident 7p centre-of-mass energies W of 134-285 GeV, has been mea- 
sured in the ZEUS detector at HERA, using an integrated luminosity of 38.4 pb~^. In- 
clusive cross sections for ep ^ 7 + X have been presented as a function of Ej' for photons 
in the pseudorapidity range —0.7 < rj'^ < 0.9, and as a function of rj'^ for photons with 
5 < Ej. < 10 GeV. The latter results have been given also for three subdivisions of the y 
range. All kinematic quantities are quoted in the laboratory frame. 

Comparisons have been made with predictions from leading-logarithm parton-shower 
Monte Carlos (PYTHIA and HERWIG), and from next-to- leading-order parton- level cal- 
culations. The models are able to describe the data well for forward (proton direction) 
photon pseudorapidities, but are low in the rear direction. None of the available variations 
of the model parameters was found to be capable of removing the discrepancy with the 
data. The disagreement is strongest in the W interval 134-170 GeV, but not seen within 
the measurement acceptance for W > 212 GeV. This result, together with the disagree- 
ments with NLO predictions seen also in recent dijet results at HERA [Q, would appear 
to indicate a need to review the present theoretical modelling of the parton structure of 
the photon. 



Acknowledgements 

It a pleasure to thank the DESY directorate and staff for their unfailing support and 
encouragement. The outstanding efforts of the HERA machine group in providing high 
luminosity for the 1996-97 running are much appreciated. We are also extremely grateful 
to L. E. Gordon, M. Krawczyk and A. Zembrzuski for helpful conversations, and for mak- 
ing available to us calculations of the NLO prompt photon cross section and a computer 
program (K&Z). 



References 

[1] ZEUS Collab., M. Derrick et al., Phys. Lett. B322 (1994) 287; 
ZEUS Collab., M. Derrick et al., Phys. Lett. B348 (1995) 665; 
ZEUS Collab., M. Derrick et al., Phys. Lett. B384 (1996) 401. 

[2] HI Collab., T. Ahmed et al., Phys. Lett. B297 (1992) 205; 
HI Collab., T. Abt et al., Phys. Lett. B314 (1993) 436; 
HI Collab., T. Ahmed et al., Nucl. Phys. B445 (1995) 195; 
HI Collab., S. Aid et al., Phys. Lett. B392 (1997) 234. 

[3] ZEUS Collab., J. Breitweg et al., Eur. Phys. J. C4 (1998) 591. 

[4] ZEUS Collab., J. Breitweg et al., Eur. Phys. J. CI (1998) 109; 

ZEUS Collab., J. Breitweg et al., DESY 99-057, to appear in Eur. Phys. J. 

[5] ZEUS Collab., J. Breitweg et al., Phys. Lett. B413 (1997) 201. 

[6] ZEUS Collab., M. Derrick et al., Phys. Lett. B297 (1992) 404. 

[7] A. Andresen et al., Nucl. Instr. Meth. A309 (1991) 101; 
A. Bernstein et al., Nucl. Instr. Meth. A338 (1993) 23; 

A. Caldwell et al., Nucl. Instr. Meth. A321 (1992) 356. 

[8] N. Harnew et al., Nucl. Instr. Meth. A279 (1989) 290; 

B. Foster et al., Nucl. Phys. B (Proc. Suppl.) 32 (1993) 181; Nucl. Instr. Meth. 
A338 (1994) 254. 

[9] ZEUS Collab., M. Derrick et al., Z. Phys. C65 (1995) 379. 

[10] R. Brun et al., CERN-DD/EE/84-1 (1987). 

[11] H.-U. Bengtsson and T. Sjostrand, Comp. Phys. Commun. 46 (1987) 43; 
T. Sj5strand, CERN-TH.6488/92. 

[12] G. Marchesini et al, Comp. Phys. Commun. 67 (1992) 465. 

[13] A. D. Martin, R. G. Roberts and W. J. Stirling, Phys. Rev. D50 (1994) 6734. 

[14] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D46 (1992) 1973. 

[15] ZEUS Collab., J. Breitweg et al., Eur. Phys. J. C6 (1999) 67; 
ZEUS Collab., M. Derrick et al., Z. Phys. C72 (1996) 399. 

[16] L. E. Gordon, Phys. Rev. D57 (1998) 235 and private communication; 

L. E. Gordon, liep-ph/9706355| and Proceedings, Photon 97, Egmond aan Zee, eds. 



A. Buijs and F. Erne (World Scientific, Singapore, 1998), 173. 



[17] L. E. Gordon and W. Vogelsang, [hep-ph/9606457| and Proceedings, DIS96, Rome, 



eds. G. D'Agostini and A. Nigro, (World Scientific, Singapore, 1997), 278. 



[18] M. Krawczyk and A. Zembrzuski, |hep-ph/9810253| and Proceedings, Photon 97, 
Egmond aan Zee, eds. A. Buijs and F. Erne (World Scientific, Singapore, 1998), 162 
and private communication. 



[19] M. Klasen, |hep-ph/9907366| . 



[20] B. L. Combridge, Nucl. Phys. B174 (1980) 243. 

[21] L. E. Gordon and J. K. Storrow, Nucl. Phys. B489 (1997) 405. 

[22] P. Aurenche et al., LPTHE preprint 92/13; 
P. Aurenche et al., Z. Phys. C64 (1994) 621. 



h/rp 


da/dE:} pb GeV"^ 


5.0- 6.0 


18.4 ± 2.1 til 


6.0 - 7.0 


9.9 ± 1.3 t\l 


7.0 - 8.0 


8.7 ± 1.1 1°:? 


8.0-9.5 


3.3 ± 0.6 til, 


9.5- 11.0 


2.2 ± 0.4 1°J 


11.0- 13.0 


1.3 ± 0.3 Igl 


13.0- 15.0 


0.3 ± 0.3 til 



Table 1: Differential cross sections for inclusive photoproduction of isolated photons 
with —0.7 < 1]'^ < 0.9, averaged over given transverse-energy intervals, for 0.2 < y < 0.9 
(134 < W < 285 GeV). The first error is statistical, the second is systematic. 



^7 


da/drj^ pb 


0.2 < y < 0.9 


0.2 <y< 0.32 


0.32 < y < 0.5 


0.5<y <0.9 




134 < W^ < 285 


134 <W < 170 


170 <W < 212 


212 < W < 285 


-0.6 


38.9 ± 5.9 tt-^7 


10.0 ± 3.3 tli 


15.2 ± 3.6 tli 


13.7 ± 3.2 tli 


-0.4 


40.1 ± 5.7 It? 


17.0 ± 3.6 tli 


13.8 ± 3.2 tli 


9.3 ± 2.8 tli 


-0.2 


27.7 ± 5.0 tli 


11.4 ± 3.1 tli 


11.7 ±2.9 1^:^ 


4.6 ± 2.1 tli 


0.0 


35.1 ± 5.5 tti 


17.7 ± 3.5 tli 


9.1 ± 3.0 t\i 


8.3 ± 2.7 tli 


0.2 


21.0 ± 4.4 tli 


9.9 ± 2.5 tli 


5.7 ± 2.4 tli 


5.4 ± 2.4 tli 


0.4 


18.7 ± 3.7 tli 


7.2 ± 2.1 tli 


7.9 ± 2.4 tli 


3.6 ± 1.7 til 


0.6 


UA ± 4.4 tli 


6.8 ± 2.7 tl:t 


5.8 ± 2.5 tli 


1.8 ± 2.2 +[!:? 


0.8 


19.5 ± 4.9 tli 


6.8 ± 2.9 t\i 


10.1 ± 3.4 tli 


2.6 ± 2.2 tli 



Table 2: Differential cross sections per unit pseudorapidity for inclusive photoproduction 
of isolated photons with 5 < E^ < 10 GeV, averaged over laboratory pseudorapidity 
intervals of ±0.1 about the given central values. The 7p centre-of-mass energy (W) 
ranges are in GeV. Results are listed for the full range of fractional incident photon 
energy y and in three subdivisions. The first error is statistical, the second is systematic. 
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Figure 1: (a) Distribution of <5Z> for prompt photon candidates in selected events, (b) 
Distribution of fmax for prompt photon candidates in selected events after cutting on <6Z>. 
Also given in both cases are fitted Monte Carlo distributions for photons, vr*^ and t] mesons with 
similar selection requirements as for the observed photon candidates. Samples with fmax > 0.75 
and fmax < 0.75 are enriched in the photon signal and in the meson background, respectively. 
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Figure 2: Differential cross section da/dEj, for prompt photons produced over —0.7 < rf < 
0.9. The inner (thick) error bars are statistical; the outer include systematic errors added in 
quadrature. The data points are plotted at the respective bin centres (see Table 1; bin centring 
corrections are negligible). Predictions are shown from PYTHIA and HERWIG at the hadron 
level (histograms), and from LG and K&Z (curves). In K&Z, the default 4-flavour NLO A-p^ 
value of 320 MeV is used. 
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Figure 3: Differential cross-section da/drj'^ for isolated prompt photons with 5 < Ej, < 10 GeV, 



for 0.2 < y < 0.9 (134 < W < 285 GeV). The inner (thick) error bars are statistical; the outer 
include systematic errors added in quadrature. Also plotted are (a) PYTHIA and HERWIG 
predictions using the GRV(LO) photon parton densities; (b) LG and K&Z NLO predictions 
using GRV(HO) and GS(HO) photon parton densities. 
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Figure 4: Differential cross-section da/drj'^, for isolated prompt photons with 5 < Ej, < 10 
GeV, compared with PYTHIA and with LG and K&Z NLO predictions, using GRV photon 
parton densities as in Fig. 3. The inner (thick) error bars are statistical; the outer include 
systematic errors added in quadrature. The plots correspond to the W ranges (a) 134-170 GeV, 
(b) 170-212 GeV, (c) 212-285 GeV. 
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